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Abstract 

We review the calculations of form factors and coupling constants in vertices with charm mesons 
in the framework of QCD sum rules. We first discuss the motivation for this work, describing 
possible applications of these form factors to heavy ion collisions and to B decays. We then present 
an introduction to the method of QCD sum rules and describe how to work with the three-point 
function. We give special attention to the procedure employed to extrapolate results obtained 
in the deep euclidean region to the poles of the particles, located in the time-like region. We 
present a table of ready-to-use parametrizations of all the form factors, which are relevant for the 
processes mentioned in the introduction. We discuss the uncertainties in our results. We also give 
the coupling constants and compare them with estimates obtained with other methods. Finally 
we apply our results to the calculation of the cross section of the reaction J/ip + ir^-D + D*. 

Prepared for publication in Prog. Part. Nucl. Phys. 
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1 Introduction 



When we write the high energy electron-proton scattering cross section using the Feynman rules of 
quantum electrodynamics (QED), we do not know how to write the interaction between the photon and 
the proton when the latter is an extended and composite object. In order to organize and parametrize 
our ignorance making use of general principles, such as charge conservation, we introduce functions of 
the involved kinematical variables, which are called structure functions. After convenient redefinitions 
these functions become form factors. The final cross sections, written in terms of form factors are 
then adjusted to data and, in this way, the form factors are measured. In a particular reference frame, 
their Fourier transform gives the spacial distribution of charged matter in the proton. This procedure 
gives a nice "picture" of the proton in space and the most accurate description of its "form" . In lower 
energy experiments, where also the four momentum transfer (g 2 ) is low, it was possible to determine 
the electromagnetic form factor (and the charge radius) of the nucleon [1]. 

QED was also applied to deep inelastic electron-proton scattering with the assumption that the 
proton is composed by pointlike constituent fermions. In this case the structure functions turned into 
the parton distribution functions and their precise measurement gave reality to quarks and gluons and 
illuminated a new world within the proton. In higher energies experiments and very large values of q 2 a 
very different picture of the nucleon emerged, in which it is made of pointlike particles, the quarks. From 
these observations one may conclude that, when probing the nucleon, nearly on-shell photons (q 2 ~ 0) 
recognize sizes whereas highly off-shell photons (q 2 0) do not [1]. This statement is supported by 
the phenomenologically very successful vector meson dominance hypothesis, according to which real 
photons are with a large probability converted to vector mesons (which are extended objects) and 
then interact with the nucleon [2]. HERA data on electron- proton reactions can be well understood 
introducing a "transverse radius of the photon" , parametrized as 

1 

where Q 2 = —q 2 and m is the mass of the vector meson considered. This empirical formula tells us 
that for Q 2 — > oo the photon is pointlike and "resolves" the nucleon target, i.e., identifies its pointlike 
constituents and does not "see" the size of the nucleon. Moreover, this formula indicates that for 
Q 2 ~ and for light mesons (like the p°) the photon has appreciable transverse radius and therefore 
also identifies the global nucleon extension. Finally, in the above formula we may have a heavy vector 
meson (J/ip or T) which will, either real or virtual, resolve the nucleon into pointlike constituents. 
This feature nicely explains why the J/ip photoproduction cross section grows pronouncedly with the 
(photon-proton) energy whereas the p cross section grows very slowly. In the former case the compact 
J/xjj interacts with the small x gluons in the protons, which have a fastly growing population. In the 
latter case, the p identifies the global and slowly growing geometrical size of the proton. 

After this glorious history as an important tool in the QED machinery, form factors were calculated 
in other contexts. In particular they were introduced in the development of effective theories to study 
interactions of charm particles. In this context, an interesting question is: which are the form factors 
of the charm mesons when probed by light particles such as photons, pions and p mesons? Apart from 
their intrinsic value as fundamental knowledge about nature, the answer to this question has immediate 
applications in hadron physics. 

Vertices with three mesons where at least two of them carry charm appear in theories of the charm 
meson interactions. This kind of theory started to become popular in the late nineties [3], during 
the analysis of the CERN-SPS data on charmonium production in heavy ion collisions. At that time 
it was believed that J/ip suppression was a good signature of quark gluon plasma formation (QGP) 
[4]. However a careful evaluation of the background was needed, in order to isolate the signal. In 
this case the background was the charmonium absorption by light mesons within the hadronic fireball 
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formed at the late stage of these collisions. Since the center of mass energy of these collisions was of 
the order of magnitude of the temperature, i.e. ~ 100 — 200 MeV, the interaction regime was clearly 
non-perturbative and the best tools were the effective Lagrangian models with charm mesons [3, 5]. 
From these Lagrangians one can derive the Feynman rules and compute scattering amplitudes. In order 
to avoid infinities one can introduce form factors in the vertices. In phenomenological applications, 
these form factors contain a parameter A which plays the role of a cut-off. In fact, even diagrams which 
give finite contributions to the cross sections must contain form factors. Otherwise the obtained cross 
sections are unacceptably large. After the introduction of form factors the results for the charmonium 
interaction cross section become quite reasonable. However these results depend too strongly on the 
choice of the cut-off parameter. In QCD sum rules (QCDSR) we can calculate these form factors from 
first principles, eliminating the freedom in the choice of parameters. 

At this point one might argue that this program is futile for at least two reasons. In first place one 
might say that one can compute cross sections such as, for example, J ftp + n — > D + D, directly from 
QCDSR [6] and it is not necessary to compute form factors. This statement is in principle correct. In 
practice however, the direct calculation of cross sections with QCDSR requires the use of the four-point 
function, which is much less precise than the three and two-point functions. It is not clear, for example, 
whether one should perform one, two or even three Borel transforms. In second place one might say that 
it is more promising to develop a chiral perturbation theory for these interactions, where one eliminates 
form factors [7]. This may indeed be the case, but for now this kind of theory is in a very preliminary 
stage. Moreover, in effective theories one needs the coupling constants. In the case of the charm three 
meson vertices these couplings (with the exception of the D*Dn coupling constant) are not measurable. 
They can be calculated in QCDSR, as a by-product of the calculation of form factors. 

The experimental study of interactions of charmed mesons with nucleons and also with light mesons 
will be one of the main topics of the scientific program of the PANDA and CBM experiments at the 
future FAIR facility at GSI [8]. 

So far we have been emphasizing the usefulness of charm form factors in the calculation of char- 
monium interaction cross sections. However, this is just one of the several applications of these form 
factors. Another context where they are needed is in heavy meson decays. Since 2003, due to the 
precise measurements of B decays performed by BELLE, BES and BABAR, charm form factors gained 
a new relevance. In B decays new particles have been observed, such as the _D s j (23 17) and the X(3872) 
[9, 10]. These particles decay very often into an intermediate two body state, which then undergoes 
final state interactions, with the exchange of one or more virtual mesons [11]. As an example of specific 
situation we may consider the decay X(3872) — >■ J/ijj + p. This decay may proceed in two steps. First 
the X decays into a D - D* intermediate state and then these two particles exchange a D* producing 
the final J ftp and p [12]. This process is illustrated in Fig. 1. In order to compute the effect of these 
interactions in the final decay rate we need the pD*D* and J/ipD*D form factors. 

Besides the decays, charm form factors may be needed to understand the structure of these newly 
observed charm states. In [13, 14, 15] it has been suggested that some of these states are dinamically 
generated resonances from the vector-vector mesonic interaction. Here we have box diagrams with 
virtual meson exchanges and vertices which need form factors. 

In this article we review our works on charm form factors [16, 17, 18, 19, 20, 21, 22, 23, 24] calculated 
with QCD sum rules [9, 25, 26, 27]. In doing so we compare our method with other approaches, 
pointing out their virtues and shortcomings. The text is organized as follows. In Section 2 we present 
a brief introduction to the QCDSR method, where we emphazise the main concepts and the strategies 
employed to calculate masses, decay constants, form factors and coupling constants. This introduction 
is intentionally abstract and "clean" . The subsequent Section 3 illustrates how to calculate numerically 
all desired quantitities. This is where the limitations of the method become clear and also where we 
show how to circunvent them. In this section we discuss the quality of numerical results and establish 
certain conditions which must be fulfilled by reliable sum rules. After having defined the method and 
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Figure 1: Final state interactions with charm meson loops. 



the quality criteria, in Section 4 we present a sample of results obtained for the form factors. We give 
special attention to one particular form factor and the respective vertex coupling constant: D*Dn. 
This vertex is special because is the only one directly accessible to experiments. Its coupling constant 
was measured and several calculations were performed, one of them with lattice QCD techniques. The 
D*Dn vertex provides the precision test for the different ways of computing the coupling constant. A 
careful analysis of this vertex resulted in the formulation of a complementary technique (the evaluation 
of meson loops in the vertex) used to improve the quality of the QCDSR calculations. In Section 5 
we show the full compilation of all the already obtained charm form factors, with the corresponding 
parametrizations, which are ready to be easily employed in phenomenological analyses of vertices with 
charm mesons. We also present a quantitative study of the uncertainties in our results. In Section 6 
we perform a comparison with the results obtained with other approaches. In Section 7 we show an 
application of the charm form factors to the study of J ftp production and absorption by light mesons. 
Finally, in Section 8 we present a summary of our results. 



QCD sum rules have been discussed in many reviews [9, 25, 26, 27] emphasizing various aspects of 
the method. QCDSR are a powerful tool to extract qualitative and quantitative information about 
hadron properties. In this approach we start with a correlation function that is constructed in terms 
of hadronic currents, which are chosen so as to have the quantum numbers of the hadrons in question. 
The basic idea of the formalism is to approach the bound state problem in QCD from short distances 
and move to larger distances, including non-perturbative effects and using some approximate procedure 
to extract hadronic masses and couplings. Although we are mostly concerned with the three-point 
correlation function, sometimes we also need two-point correlation functions. For this reason and also 
in order to illustrate some basic procedures in QCDSR we devote the next subsections to the discussion 
of two-point correlators. 

2.1 The two-point correlation functions 

A generic two-point correlation function is usually written as: 



2 QCD Sum Rules 




(2) 



where j(x) is a current with the quantum numbers of the hadron which we want to study. 
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In QCDSR we assume that the correlation functions may be written at both the quark and the 
hadron levels. Identifying the hadronic representation with the corresponding representation in terms 
of quarks and gluons, we obtain the sum rule, from which we can determine hadron properties. In 
the QCD side we proceed with the calculation of the correlation function using the operator product 
expansion (OPE), which is formulated with Wilson coefficients and local operators in terms of the non- 
perturbative structure of the QCD vacuum. In order to apply this method to the correlation function 
(2), we have to expand the product of two currents in a series of local operators: 

n OP V) = i I d 4 xe^(0\T[j(x)j\0)\0) =J2C n (q 2 )(0\d n \0}, (3) 

J n 

where the set {O n } includes all local gauge invariant operators expressible in terms of light quark and 
gluon fields. By construction the coefficients C n (Q 2 ) (Q 2 = —q 2 ) in Eq. (3) include only short-distance 
effects. Therefore, they can be evaluated perturbatively. Non-perturbative long-distance effects are 
contained only in the local operators. In the OPE expasion, the operators are ordered according to 
their dimension n. The lowest- dimension operator with n = is the unit operator associated with 
the perturbative contribution: C (Q 2 ) = W er (Q 2 ), O = 1. The QCD vacuum fields are represented 
in (3) in the form of the so called vacuum condensates. The lowest dimension condensates are the 
quark condensate of dimension three: (0|O3 10) = (qq), and the gluon condensate of dimension four: 
(0|O4|0) = (g 2 G 2 ). The contributions of higher dimension condensates are suppressed by large powers of 
Aq CD /Q 2 , where Aqcd is the typical long-distance scale. Therefore, even at intermediate Q 2 ~ 1 GeV 2 , 
the expansion in Eq. (3) can be safely truncated after a few terms. 



2.2 The spectral density 

The hadronic (or "phenomenological" ) representation of the correlation function in Eq. (2) can be 
written as a dispersion relation: 

W hen (q 2 ) = - fds - P(S) , + • • • , (4) 
J q l — s + le 

The dots in the above equation represent subtraction terms. The spectral density is given by the 
imaginary part of the correlation function: 

p(s) = -Im[U(s)} . (5) 

7T 

The evaluation of p(s) is simpler than the evaluation of the correlation function itself, and the knowledge 
of p(s) allows one to recover the whole function IT(g 2 ) through the integral in Eq. (4). 

The current j (j^) appearing in (2) and (3) is an operator that annihilates (creates) all hadronic 
states that have the same quantum numbers as j. Consequently, IT(g 2 ) contains information not only 
about the low mass hadron of interest, but also about all excited states with the same quantum numbers. 
When comparing a set of hadrons with the same quantum numbers, the lowest resonance is often fairly 
narrow, whereas higher mass states are broader. We can therefore parametrize the spectral density as a 
single sharp pole, representing the lowest resonance of mass m, plus a smooth continuum, representing 
higher mass states: 

p(s) = X 2 5(s-m 2 )+p cont (s), (6) 

In the above Equation the parameter A represents the coupling of the current to the low mass hadron 
H: 

(0\j\H) = \. (7) 
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For simplicity, we often assume that the continuum contribution to the spectral density, p cont (s) in 
Eq. (6), vanishes bellow a certain continuum threshold Sq. In order to keep the number of parameters 
as small as possible, we assume that above s the spectral density is given by the result obtained with 
the OPE. This idea is implemented by the Ansatz proposed in [28] : 

p^(s)=p OPE (s)Q(s-s ). (8) 

This is called quark-hadron duality. 
2.3 The mass sum rule 

The sum rule is obtained from the matching of the two descriptions of the correlator: 

U phen {Q 2 ) = U OPE {Q 2 ) . (9) 

However, such a matching is not yet practical. The phenomenological description is significantly dom- 
inated by the lowest pole only for sufficiently small Q 2 , or even better, timelike q 2 near the pole. On 
the other hand, the OPE side is only valid at a sufficiently large spacelike Q 2 . In order to improve the 
overlap between the two sides of the sum rule, we apply the Borel transform: 

/ 2\n+l / j \ n 

£ M2 [I% 2 )]= Jim LJLL * n(g 2 ). (10) 

-q 2 /n=M 2 

It is interesting to notice that: 

B M 2 [g 2 ™] =0, n > . (11) 

This means that all the subtractions terms in Eq. (4) are eliminated by the Borel transform. Another 
important feature of the Borel tranform is the fact that: 



1 



(m 2 — q 2 ) n 



Y e ~m?/M 2 

n>0. (12) 



{n-l)\ (M 2 )™- 1 ' 



Therefore, the Borel transform exponentially suppresses the contribution from excited resonances and 
continuum states in the phenomenological side. In the OPE side the Borel transform suppresses the con- 
tribution from higher dimension condensates by a factorial term, improving the OPE convergence. After 
making a Borel transform on both sides of the sum rule, and transferring the continuum contribution 
to the OPE side, the sum rule can be written as 

x 2 e -m 2 /M 2 = f S ° ds e -s/M 2 p OPE {s) (13) 

J s min 

The mass of the low-lying state, m, can be determined by taking the derivative of Eq. (13) with respect 
to 1/M 2 , and dividing the result by Eq. (13). This gives: 

f s ° ds e~ s,M ' 2 s o OPE (s) 

J S S0 ds e~ s / M2 p OPE (s) • 1 ' 

Since in the evaluation of both sides of the sum rule we have to make approximations, the value extracted 
from Eq. (13) will be a function of M 2 . The Borel window is defined as the range of values of M 2 where 
the two sides of the sum rule have a good overlap and, therefore, information on the lowest resonance 
can be extracted. In general the Borel window is determined by imposing two different criteria: the 
minimum value of the Borel mass is fixed by requiring the convergence of the OPE and the maximum 
value of the Borel mass is determined by imposing the condition that the pole contribution must be 
bigger than the continuum contribution 
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2.4 The three-point function 

The three-point function associated with a generic vertex of three mesons M 1; M 2 and M 3 is given by 

r(p,p') = J d A xd A y e^*e^^(0\T{j s (x)jl(y)A(0)}\Q) , (15) 

where the current ji represents states with the quantum numbers of the meson %. As in the previous 
Section, the correlation function is evaluated in two ways. In the first one, we consider that the 
currents are composed by quarks and write them in terms of their flavor and color content with the 
correct quantum numbers. This is the QCD description of the correlator and it gives rise to the QCD 
side (or OPE side) of the sum rule. In the second way, we write the correlation function in terms of 
matrix elements of hadronic states which can be extracted from experiment, or calculated with lattice 
QCD or estimated with effective Lagrangians. In this last approach we never talk about quarks and use 
all the available experimental information concerning the masses and decay properties of the relevant 
mesons. This is the hadronic description of the correlator and is called the phenomenologial side of the 
sum rule. After studying both sides separately, we identify one description with the other and write 
the sum rule. 

2.5 The OPE side 

We start with the current, which has the general form: 

Ji = ?Ag, (16) 

with A = 1,7^,75,7^75 for a scalar, vector, pseudo-scalar and axial-vector meson, respectively, q is the 
quark spinor field. Since the currents may carry Lorentz indices, so will the vertex function T. When 
we insert the three currents into Eq. (15) we get the vacuum expectation value of the T product of 
six quark fields multiplied by Dirac matrices in different points. In this expression we apply Wick's 
theorem obtaining a series of terms, each of which being the product of contractions (propagators) 
times a vacuum expectation value of normal ordered operators taken in the QCD vacuum. These latter, 
in the local approximation, are the QCD condensates. This series is precisely (3) and its first terms 
are diagramatically depicted in Fig. 2. The first and leading one is the perturbative term in Fig. 2a. 
Figs. 2b and 2c show examples of terms with the quark condensate and with the gluon condensate 
respectively. 




Figure 2: Diagrammatic representation of the first terms of the OPE series for a three-point funtion. 
(a) The perturbative term, (b) The quark condensate, (c) The gluon condensate. 

The evaluation of (15) in lowest order leads to a loop diagram with three quark propagators. The 
first corrections to the simple bubble (called the perturbative contribution) come from diagrams where 
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a quark line "dives" into the vacuum and emerges from it. This contribution to the correlators appears 
multiplied by the quark condensate (qq) . Higher order corrections represent the interactions of the quark 
lines exchanging gluons among themselves (a s corrections) and with the vacuum (gluon condensates and 
quark-gluon mixed condensates). An equivalent description is obtained performing first the contractions 
in (15), obtaining three quark propagators and then performing an operator product expansion (OPE) 
of these propagators. For our purposes it is enough to consider only the first few terms of this expansion: 

S ab (x) = (0|T[g a (x)g fe (0)]|0) = ^-J - £0 - <*Z> + ^H<^ 

T^GU^ + <^t) - 7^X b 9sG±^ H-x 2 ). (17) 



32n 2 x 2 ab ^~^ r ~ ~ r > 32tt 2 aW ^ 
In the above expression x is the four-vector which defines the separation between the two quark fields, 
m q is the quark mass, are the SU(3) generators, G^ v is the gluon field tensor and a^ v is the 
anticommutator of the Dirac matrices. In (17) the first two terms are related to Fig. 2a, the next two 
terms are related to Fig. 2b and the last two are related to Fig. 2c. After the contractions and the 
OPE expansion we arrive at an expression for T which has the generic form: 

roP£ = E^ OP V^V)^, (is) 

3 

where F° PE (p 2 , p' 2 , q 2 ) are invariant functions of the momenta and Lj are the structures, i.e., products 
of Dirac matrices, the metric tensor and the four momenta, carrying Lorentz indices. For each one 
of the invariant amplitudes appearing in Eq.(18), we can write a double dispersion relation over the 
virtualities p 2 and p' 2 : 

F OPE = - J- f ds r du P° PE ^ U '^ , (19) 

47T 2 Js mm Ju mm (S - p 2 )(u - p' 2 ) ' 

where pf PE (s, u, q 2 ) is the double discontinuity of the amplitude Ff PE (p 2 ,p' 2 , q 2 ) and can be calculated 
using the Cutkosky's rules. We can work with any structure appearing in Eq.(18), but we must choose 
those which have less ambiguities in the QCD sum rules approach, which means among other things, a 
weak influence from the higher dimension condensates. The invariant amplitudes and thus the double 
discontinuities p, receive contributions from all terms in the OPE. The first one (and dominating) of 
these contributions must come from the perturbative term. 

2.6 The phenomenological side 

For the phenomenological side we formally repeat the steps mentioned above and find that, as before, 
the correlation function can be written as a sum of contributions with different tensor structures, i.e. : 

T vhen = Y.Ff en {v\v\f)^ (20) 

3 

where Fj hen (p 2 ,p' 2 , q 2 ) are invariant functions of the momenta and Lj are the structures. For each one 
of the invariant amplitudes appearing in Eq.(20), we can write a double dispersion relation over the 
virtualities p 2 and p' 2 , holding Q 2 = —q 2 fixed: 

1 r°° r°° o phen (s u O 2 ) 
Ff h ™ = ~ ds du 2 \, L (21) 

47T 2 Js mm Ju mm [S - p Z )(u - p' ) 

where p phen (s,u,Q 2 ) is the double discontinuity of the amplitude F[ hen (p 2 , p' 2 , Q 2 ). The function p can 
be generically written as: 

p? hen (s,u,Q 2 ) = a5(s — m 2 ) 5(u — to 2 ,) + b5(s — to 2 ) 9{u — uq) + 

c5(u-m 2 2 )9(s-s ) + p™ nt (s,u,Q 2 ) 9(s - s ) 9(u - u ), (22) 
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where sq and uq are the continuum threshold parameters. In the above equation the first term represents 
a kinematical situation where the mesons M x and M 2 (with masses mi and m 2 respectively) are in the 
ground state and M 3 is off-shell, having an arbitrary Euclidean four momentum squared Q 2 . The second 
term represents a situation where the meson Mi is on-shell but the groundstate of the meson M 2 is 
absent in the vertex, which contains only its excitations starting at uq. The third term represents the 
inverse situation, where M 2 is on-shell and it is accompanied in the vertex by the excitations of Mi 
starting at sq. Finally, the last term represents the excitations of Mi and M 2 , which start at sq and 
uq respectively. After a double Borel transform in the variables p 2 and p' 2 the second and third terms 
are exponentially suppressed. One can then safely assume quark-hadron duality and parametrize the 
continuum by the double discontinuity of the theoretical part: pf ont (s,u,Q 2 ) = p° PE (s,u,Q 2 ). This 
point is thoughtfully discussed in [30]. Inserting the above expression into (21) (and neglecting the 
suppressed terms) we have: 

- A -"" - h i: d ° /; *• <23) 

where A? hen refers to the ground state (or pole) contribution in the structure i. In order to calculate 
hf Len we go back to (15) and insert complete sets of hadronic states. After some algebra [31] we obtain: 

j^phen _ !mJmJm z (M 1 M 2 M 3 ) i 
1 (p 2 — m 2 )(p' 2 — m\){q 2 — m 2 ) 

The meson decay constants appearing in the above equation are defined by the following matrix ele- 
ments: 



(0\%\V)=m v f v e„ 


(25) 


(0\j£\A) =m A f A € li , 


(26) 


(o\MP) = —fp, 


(27) 


m q 


(0\f A \P)=ifpp», 


(28) 



and 

for vector, axial, pseudoscalar and axial currents respectively. In the above expressions e M and p M are 
the relevant polarization vector and four momentum of the involved mesons resepectively. Eq. (28) 
refers to the case where a pseudoscalar meson is represented by an axial current, and m q in Eq. (27) is 
the heaviest quark in the pseudoscalar meson P. 

Since the higher states have been considered in the second term of (23), the amplitude (M 1 M 2 M 3 ) 
appearing in (24) refers only to the ground states of the mesons Mi, M 2 and M 3 and can be calculated 
with the help of an effective Lagrangian of the type: 

£m 1 m 2 m 3 = 9m 1 m 2 M:, (v AMi AM 2 M 3 + hcj , (29) 

where r\ = i or e^ ua/ 3, A = 1 or <9 M , he stands for Hermitian conjugates and gM 1 M 2 M 3 is the coupling 
constant. The choices implied in r\ and A depend on the specific combination of mesons. 

2.7 The sum rule 

The sum rule is obtained from the identity: 

F l phen (p,p',q) = F l OPE (p,p',q), (30) 
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that gives 



/mi /m 2 /m 3 ( Mi M 2 M 3 ) j 1 r° d „ T° iu pP Pg (a,«,Q 2 ) (31) 



(p 2 - mf)(p' - m|)(g 2 - m§) 47r 2 A min A min (s - p 2 )(u - p' 



where we have used (19), (23) and (24) and have transfered the dispersion integral in the phenomeno- 
logical side to the OPE side. Notice that the integration limits are now finite. In order to improve the 
matching between the two sides of the sum rule and also to suppress the pole-continuum transitions 
[30] we perform a double Borel transform (10) in the variables P 2 = —p 2 — > M 2 and P' 2 = —p' 2 — > M' 2 , 
on both sides: 

fMjMMMMM^ e „ m?/M2 e _ m|/M , 2 = 1 fso ^ ruo dup OPE {s ^ Q 2 )e -s/M> e -u/M^ ^ 

[Q +m 3 ) Arc Js min Ju min 

From (29) it is clear that the above equation can be solved for the coupling gM 1 M 2 M 3 - Since the squared 
momenta of the mesons Mi and M 2 (p 2 and p' 2 respectively) have been replaced by the Borel masses M 2 
and M' 2 and then fixed, the coupling g will be a function only of the remaining Euclidean momentum 
Q 2 , i.e., g ( M 1 3 M 2M3 = 5 , m^m 2 m 3 ( ( 3 2 ) an d this is what we call a form factor. The superscripts in parenthesis 
indicate the meson M 3 is off-shell. At the point Q 2 = —m 2 the meson M 3 is on-shell and the form 
factor becomes the coupling constant. Of course one can not use Q 2 = —m 2 in Eq. (32) since the sum 
rule is only valid for Q 2 > 0. Therefore, to obtain the coupling constant we will need to use some 
extrapolation procedure, that will be discussed in Sec. 4.1. 



2.8 Effective Lagrangians 

Since the pioneering work of Matynian and Miiller [3] , there has been an intense discussion concerning 
the details and properties of the effective Lagrangians which describe the interactions among charm 
mesons. Here we follow Refs. [32, 33]. We write down the SU(4) chiral Lagrangian in terms of the 
mesonic fields: 

C = Tr (d^PP) - ^Tr {F%F^) , (33) 

where = dpV v — d v V^ and P and V denote respectively the properly normalized 4x4 pseudoscalar 
and vector mesons matrices in SU(4) given by: 



V = 



1 

72 



1 

72 



t ^ + 



V 

/ Pi 



V 



-2- 4 


Vc 




7T + 




K + 




D° \ 








V2 






K° 




D~ 




K~ 






K° 




fiv + : 


Vc 

/12 


D- 


■ 


D° 






D+ 




Df 


3?2c 

Vu / 




juL + 


il± 


P" 


P + 




K *+ 




D*° 




p~ 




+ 7/6 + 


m 


K*° 




D*- 










K*° 


~\ 


fl"' + 


JJ± 
Vl2 


D*~ 




D*° 






D *+ 




D*+ 


3J/tp 


) 



(34) 



(35) 



In order to obtain the interaction Lagrangian we introduce the following minimal substitutions: 



d,P -4 d,P- l j[V„P] 



(36) 
(37) 
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After using the Hermiticity of P and V the resulting Lagrangian reduces to: 

C = C + tgTT(d^P[P,V,])- 9 jTr([P,V,] 2 ) 

+igTr ( W [V„ V v \) + ^Tr ([V„ V u f) . (38) 

The above Lagrangian accounts for vertices of the form PPV and VVV, where P and V denote 
pseudoscalar and vector mesons respectively. In order to include the vertices of the form PVV, it is 
necessary to use the following anomalous three-particle Lagrangian 



£L = -0rf ua " TT • ( 39 ) 

Equations (38) and (39) are used to write the relevant meson effective Lagrangians. For example, in 
the PPV (VVV) case, we choose two (none) mesons in (34) and one (three) meson in (35) and set 
all the other elements of these matrices to zero. We then substitute the resulting Equations (34) and 
(35) into (38) and obtain the interaction Lagrangian for one specific three meson interaction vertex. 
The coupling constants appearing in these Lagrangians are charge specific (for example: gD+D-j/ip) and 
are functions of the universal SU(4) couplings g and g a , being thus connected to each other. They are 
listed in Table 1 and their connection with g and g a is shown. Using a compact notation, the relevant 
Lagrangians are: 

C*d*d = tg^DD;(Dd^-d^Dn), (40) 
C^d*d = gip D *De^d^ u (d a D;D + Dd a D*p), (41) 
C^dd = ig+DDV (?»DD - DdpD) , (42) 

CvD'D* = -gvD-D-^dnDydaD'p, (43) 

C^d*d* = ig^*D* [V (dpD^Dt - D^d^D*) , 

+ (d fl ip u D* u - ^Dl) D*» + D*^ (rd,D* u - d^D*")] , (44) 

C p dd = -ig pDD ff (d^DD - Dd^D) , (45) 
C P D*D* = ig P D*D* y (dpD* v D* v - D* v dpD*) , 

+ (d^D*" - p v d„Dl) D*» + D*» (p u d,D* u - d, Pv D*») ] , (46) 

CpD*D = g P D*D^ aP d,p v (d a Dp + Dd a D* p ) (47) 

For future applications we present below the following Lagrangians containing quartic vertices: 

C^DDn = ig^DDn^^Dd^dpD, (48) 

Ci,D*D* = -gwDvV (DttD% + D;tiD) , (49) 
C^d*. = tg^D*D*,e^^D\d a nD*p + ih^ D * D * n e^d^D* a nD*p, (50) 

where it is understood that, when doing practical calculations, charges must be specified in the above 
expressions. As it will be seen in the next section, after having calculated all specific coupling constants, 
we can, using some appropriate convention, write them all in terms of generic coupling constants 

From the above Lagrangians we can derive the Feynman rule for the vertex (the amplitude (M1M2M3)) 
which, as it can be seen, depends on the relevant polarization vector, momenta and on the coupling 
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gM 1 M 2 M 3 - The latter is the unknown, which will be determined by solving the sum rule. For example, 
for the D*Dn vertex we get: 

(D*(p)\n(q)D(p-q)) = sWA^) , (51) 

where the momentum assignment is specified in the brackets and e M is the polarization vector of the 
D*. 

2.9 Couplings 



Coupling 


Normal Vertices 


i//(4V3) 


ujD°D°, uD+D-,uD*°D*°, ujD*+D*- 




t]D*°D*°,7]D ± D*^ 


<?/4 


n D°D* ,n D ± D*T 




p°D°D°, p°D + D~, p°D*°D*°, p°D* + D*- 




ir + D°D*-, ir-D°D*+, ir+D-D* , n-D+D* 




p+D°D-,p-D°D + , p+D*°D*-,p-D*°D*+ 


9/VQ 


?7 cj D j D* ,77 cj D ± j D*t 




J/ipD°D°, J/i/jD+D-, J/ipD*°D*° : J/tpD^D*^ 


Anomalous Vertices 


a /(Ay/6) 


Vcp°P°,Vcp + P ,1 c ww 




J/im°p°, J/ijm^p*, J/ip7]u 


«/(4 v / 3) 


riutu, r]D* D* , r]D* + D*~ ,ujD°D*°, uD°D*°, uD ± D*^ 


«/(2 v / 6) 


r] c D*°D*°, r] c D* + D*~,r] c J/^J/tlj 




J/ijD°D*°, J/^D°D*°, J/i)D ± D*T 


a/4 


n 0D*°D*°, n D* + D*-,p D°D* , p° D° D*° , p° D ± D* T 


a/(2 v / 2) 


p +D°D*-,p-D°D* + , p+D-D*°, p-D + D*° 



Table 1: Coupling constants for the charmed three- meson vertices within the SU(4) scheme, a = 
g 2 a N c /(16rr 2 F n ). 

Using the formulas developed in Section 2.8 we can write all the specific coupling constants in terms 
of the universal couplings g and g a . The SU(4) scheme is shown in Table 1, with the help of which 
we can relate all the coupling constants of charged states among themselves. Therefore, using one 
known coupling constant, such as, for example the experimentally measured gD*+D°ir-i we can infer 
all the others. The couplings determined in this way are already enough to be used as input in the 
Lagrangians (40) - (47), which are then used in calculations of phenomenological interest. However, in 
the literature, one often finds a "generic" coupling constant, i.e. without any charge specification, such 
as, for example gD*Dw The choice of the generic go'Dir (among the the possible charged state couplings 
gD*+D°TT-i gD*°D°n°i •••etc) is a matter of convention, as mentioned explicitly in [3]. For the purpose of 
comparison we define the generic coupling as follows. If there is an isovector meson in the vertex, the 
generic coupling is the one with a neutral isovector. For example: 

gD*Dw = fi , D* ± J DT 7 r = gD*°D°TT° 

From Table 1 the generic coupling can be obtained from the states with the charged isovector: 

1 1 

gD*Dir — —7EgD*-D°n+ ~ g D* + D"tt- 
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If there is no isovector in the vertex, all the states have the same coupling and there is no ambiguity. 
Then, for example: 

gj/iPDD = gj/TpD°D° — 9J/4>D+D- 

All the calculations discussed throughout this review were performed for charged isovector currents. 

In the next section we shall discuss in more detail the criteria which must be satisfied by the sum 
rules, in order to be reliable. 



3 Evaluation of the sum rules 
3.1 Numerical inputs 

In the following subsections we will present numerical results. In the quantitative aspect, QCDSR is not 
like a model which contains free parameters to be adjusted by fitting experimental data. The inputs 
for numerical evaluations are the following: i) the vacuum matrix elements of composite operators 
involving quarks and gluons which appear in the operator product expansion (3). These numbers, 
known as condensates, contain all the non-perturbative component of the approach. They could, in 
principle, be calculated in lattice QCD. In practice they are estimated phenomenologically. They are 
universal and, once adjusted to fit, for example, the mass of a particle, they must have always that 
same value. They are the analogue for spectroscopy of the parton distribution functions in deep inelastic 
scattering [1]; ii) quark masses, which are extracted from many different phenomenological analyses; in) 
the threshold parameter s is the energy (squared) which characterizes the beginning of the continuum, 
as shown in (8). Typically the quantity ^/so — m (where m is the mass of the ground state particle) 
is the energy needed to excite the particle to its first excited state with the same quantum numbers. 
This number is not well known, but should lie between 0.3 and 0.8 GeV. If larger deviations from this 
interval are needed, the calculation becomes suspicious. All in all, in QCDSR we do not have much 
freedom for choosing numbers. In the calculations discussed below we will use the numerical inputs 
shown in Table 2. The masses and decay constants are taken from the literature [34]. They might also 
be calculated with QCDSR as discussed in previous subsections. 



Parameter Value 

m c ( GeV) 1.27 ± 0.1 

m D (GeV) 1.86 

m D *(GeV) 2.01 

m p (GeV) 0.775 

ra^GeV) 3.1 

fj/^GeV) 0.405 ±0.015 

/ D *(GeV) 0.24 ±0.02 

/o(GeV) 0.18 ±0.02 

f p (GeV) 0.160 ±0.005 

/ w (GeV) 0.131 ±0.001 

{qq)(GeVf (-0.23 ± 0.01) 3 

{g 2 G 2 )(GeV) 4 0.88 ±0.3 

Table 2: Parameters used in the calculation with their errors. 



3.2 Borel stability 

As it was mentioned above we perform two Borel transforms introducing the two Borel masses M 2 
and M . The sum rules, expressing the interesting quantities as a function of the Borel parameters, 
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must be as much independent of these parameters as possible. An extensive check of this dependence 
has been carried out in [16, 17, 18, 19, 20, 21, 22, 23, 24]. For the sake of illustration we present 
below some numerical analysis taking the vertex J/ipD*D as working example. For the continuum 
thresholds we have used so = (m^ + A s ) 2 and uo = (ttid(d*) + A u ) 2 for the sum rule where D(D*) is 
off-shell, and so = (mj). + A s ) 2 and uq = (m^ + A u ) 2 for the sum rule where J/ip is off-shell. We 
take A s = A u = (0.5 ± 0.1) GeV. We first discuss the J/ipDD* form factor with an off-shell D* meson. 




Figure 3: M 2 and M' 2 dependence of gfoo^Q 2 = 0). 

Fixing Q 2 = —q 2 = and A s = A u = 0.5 GeV we show in Fig. 3 the Borel dependence of the form 
factor g^)£)*(0)- We see that we get a very good stability for the form factor as a function of the two 
independent Borel parameters in the considered Borel regions. 

Since the dependence of the form factor on M 2 and M' 2 is weak, there is no need to continue our 
analysis with these two independent variables. We can relate one to another and use only one variable. 
One commonly used relation is: 

M' 2 _m 2 M 

M*~^ (52) 

where Mi and M2 are the mesons with P 2 = —p 2 — > M 2 and P' 2 = —p' 2 — > M' 2 respectively. If the 
meson Mi is light and the meson M2 is heavy, the above equation is replaced by: 

~M 2 ~ {m 2 M2 -mlY (53) 

Other relations between M 2 and M' 22 are possible, such as: 

M' 2 = aM 2 + b. (54) 

In [35] a comparative analysis was performed leading to the conclusion that all these forms are accept- 
able. In what follows we will discuss results obtained with (52) and with (53). We show, in Fig. 4, the 
behavior of the form factors gj/^DD* (Q 2 ) at Q 2 = 2 GeV 2 as a function of the Borel mass M'. The solid 
line gives g[f/J DD * at a fixed ratio M' 2 /M 2 = m 2 D /m 2 b . The dashed line gives g^/^cc* at a fixed ratio 
M' 2 /M 2 = mf^/m^, and the dotted line gives gf/^D* at a & xe d ra ti° M /2 /M 2 = m 2 D /m 2 D *. We can 
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Figure 4: M' 2 dependence of gji^ DD * (solid line), g^jj^ DD * (dashed line) and gj/^, D * (dotted line) for 
Q 2 = 2.0 GeV 2 . 

see that the QCDSR results for g^J DD * and gu^ DD * are very stable in the interval 3 < M' 2 < 6 GeV 2 . 
In the case of g^jfff) D * the stability is not as good as for the other form factors, but it is still acceptable. 

3.3 OPE convergence 

As mentioned in the previous sections the operator product expansion should always be convergent. 
However sometimes this convergence is not so fast. In this section we present some examples, encoun- 
tered in our calculations, which illustrate extreme cases ranging from a very fast convergence, where 
the second term is already negligible, to a slower convergence. Fortunately, in these latter cases, the 
obtention of a very weak dependence of the form factor on the Borel mass (i.e., the ocurrence of a 
stability "plateau") indicates that the sum rule is still enough convergent. 

In Fig. 5 we show the perturbative (dashed line) and the gluon condensate (dotted line) contributions 
to the form factor gj/^DD*(Q 2 ) at Q 2 = 2 GeV 2 as a function of the Borel mass M' 2 at a fixed ratio 
We see that the gluon condensate contribution is negligible, when compared with 
the perturbative contribution. The same kind of stability is obtained for other values of Q 2 and for the 
other two form factors [22, 36]. In Fig. 6 we show another example of OPE behavior, now in the D*D7r 
vertex. In this Figure we study the D*Dir vertex, choosing the pion to be off-shell with Q 2 = 1 GeV 2 
and setting M' 2 = M 2 . The solid line shows the sum of the first two relevant terms of the OPE, the 
perturbative term and the gluon condensate. It is reassuring to observe that the former is much larger 
than the latter and also that their sum is very flat over a wide range of values of M 2 . Although here 
the dominance of the perturbative term over the gluon condensate is less pronounced, we observe an 
interesting feature, namely, that the inclusion of the gluon condensate contributes to the Borel stability 
of the sum rule [16, 19]. In Figs. 7 and 8 we show a comparison between the perturbative term and 
quark condensate term. They illustrate how small (Fig. 7) and how large (Fig. 8) the effects of the 
quark condensate can be. 
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Figure 5: M' 2 dependence of perturbative contribution (dashed line) and the gluon condensate contri- 
bution (dotted line) to g^f/^ DD * at Q 2 = 2.0 GeV 2 . The solid line gives the final result for the form 
factor. 



3.4 Pole versus continuum 

In the study of the two-point function a crucial assumption is the form of the spectral density, Eq. (6). 
Moreover, it is also crucial that the dispersion integral be dominated by the pole contribution and not 
by the continuum. In the case of the two-point function, if this condition is not satisfied we can not 
compute the mass of the particle. A similar condition must hold for the three-point function. If it is 
not dominated by the pole contribution, this could be an indication that there is more than one off-shell 
particle in the vertex. In the case of the two-point function the pole dominance can be tested with the 
help of (4). The pole and continuum contributions are given by: 



Ti{q 2 ) po ie 



so P ° PE (S) 

as ■ 



q 2 — s + ie 



'cont 



ds 



p OPE {s) 



(55) 



(56) 



'so <r - s + ie 

We can perform a Borel transform in the above expressions and plot the relative contributions, which 
are given by: 

U(M 2 ) pole 



Pole 



Continuum 



Tl{M 2 ) pole + n(M 2 ) 

cont 

n(M 2 ) cont 



n(M 2 ) poie + IT(M 2 



(57) 



(5? 



'cont 



A similar analysis can be carried out for the three-point function. Starting with (19) and choosing one 
of the % tensor structures we can write: 



pOPE 
pole 



7 OPE 
cont 



4vr 2 
1 



p OPE (s,u,Q 2 ) 



i ./.so ds L du ( s - p*)(u - p> 2 y 



(59) 
(60) 



17 



II 

o 



b2 




Figure 6: D*D7r form factor as a function of the Borel mass M 2 for an off-shell pion. Dotted, dashed 
and solid lines represent the contribution of the gluon condensate, of the perturbative term and their 
sum (total) respectively. 



and then, performing a double Borel transform and using (52) or (53), we write the two relative 
contributions as: 

Pole = F OPE (M\ ole 

ForE( M 2 )pole + F oPE iM 2 )amt > K ) 

F OPE (M 2 ) 

Co,lti,mum = ?^VTf%U' (62) 

In Figs. 9 and 10, we show the pole-continuum analysis for the D*D*p vertex [23]. In Fig. 9 the 
off-shell particle is the p, with virtuality Q 2 = 1 GeV 2 . In Fig. 10 the off-shell particle is the D* with 
Q 2 = 1 GeV 2 . As it can be seen, for masses higher than ~ 2.5 GeV 2 the sum rule is dominated by the 
continuum. This condition sets an upper limit for the Borel mass. A lower limit comes from imposing 
the OPE convergence, as it can be inferred from Fig. 6. This interval is called the "Borel window" and 
it does not always exist. Changing the value of Q 2 may help or hinder the existence of the Borel window. 
In the next sections we will present QCDSR results for the form factors. We will show curves of g{Q 2 ) 
as a function of Q 2 for a fixed M 2 . The points where we have results are restricted to a certain region 
(sometimes relatively narrow) of the Q 2 axis. One of the reasons for this limited range of applicability 
is the requirement of the pole dominance. 



3.5 Continuum threshold effects 

The results showed in Figs. 4, 5, 6, 9 and 10 were obtained using A s = A u = 0.5 GeV. In Fig. 11 we use 
the form factor gjJ^ DD *{Q 2 ) to illustrate the dependence of our results on the continuum thresholds. In 
the Figure there are three sets of three curves. They show fits of the QCD sum rule results. The lower, 
intermediate and upper sets were obtained using A u = 0.4 GeV, 0.5 GeV and 0.6 GeV respectively. 
The dashed, solid and dotted lines in each set were obtained using A s = 0.4 GeV, 0.5 GeV and 0.6 GeV 
respectively. In this case we can see that the dispersion in the region < Q 2 < 4.5 GeV 2 , where we have 
the QCDSR points, does not lead to a bigger dispersion at Q 2 = —m 2 D , where the coupling constant is 
extracted. As it will be discussed in the next subsection, the extrapolation procedure used here does 
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Figure 7: D*D*p form factor as a function of the Borel mass M 2 for an off-shell D*. Dotted, dashed 
and solid lines represent the contribution of the quark condensate, of the perturbative term and their 
sum (total) respectively. 

not necessarily amplify the uncertainties. In fact, in some cases the uncertainties can be even damped 
as we move to the time-like region. An example of this damping is shown in Fig. 12, where we plot 
fits of the QCDSR results for the form factor of the D*Dp vertex. The two sets of lines correspond 
to a p off-shell (steeper lines) and to a D off-shell. Each line corresponds to a different choice of the 
continuum threshold parameters. As it can be seen, this is one of the major sources of uncertainties. 

3.6 Choice of the structure 

As it was discussed in the previous sections, the explicit evaluation of the correlation functions both 
in the OPE and in the phenomenological side leads to expressions written in terms of several tensor 
structures. We can write a sum rule identifying the coefficients of each structure and hence we have 
as many sum rules as structures. In principle all the sum rules are equivalent and should yield the 
same final results. In practice however, the truncation of the OPE affects different structures (and the 
corresponding sum rules) in different ways. Consequently some structures lead to sum rules which are 
more stable. In the simplest cases, such as in the D*Dp vertex, we have only one structure and, using 
the notation of (18) and (20), the correlator (both in the phenomenological and the OPE descriptions) 
is written as: 

r^(p,p') = F(p 2 ,p' 2 , q 2 ) e a ^p a P ,P , (63) 

and an analogous expression for the OPE side. In the most complicated cases, as in the pD*D* or 
J/ijjD*D* vertices, the number of structures is fourteen and the correlators (both in the phenomeno- 
logical and the OPE descriptions) have the following tensor decomposition: 

r»va(p,p') = F x {p 2 ,p' 2 , q 2 )g^Pa + F 2 (p 2 ,p' 2 , q 2 )g m Pu + F 3 (p 2 ,p' 2 , q 2 )g ua p^ 
+F 4 (p 2 ,p' 2 , q 2 )g^p' a + F 5 (p 2 ,p' 2 , q 2 )g m P v + F 6 (p 2 ,p' 2 , q 2 )g ua p'^ 
+F 7 (p 2 ,p' 2 , q 2 )p^p u p a + F 8 (p 2 ,p' 2 , q^p'^Pa + F 9 (p 2 ,p' 2 , q 2 )p^p' v p a 
+F 10 {p 2 ,p' 2 , q 2 )p^PuP a + F n (p 2 ,p' 2 } q 2 )pjp' u p a + F 12 (p 2 ,p 2 , q 2 )p^ v p' a 
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Figure 8: D*D*ir form factor as a function of the Borel mass M 2 for an off-shell D*. Dashed, dotted 
and solid lines represent the contribution of the quark condensate, of the perturbative term and their 
sum (total) respectively. 

+f 13 (p 2 , v '\ q 2 )pM + Fu( P 2 , v '\ q 2 W„p' a - (64) 

The above equation is the most general expression that can be written with three Lorentz indices //, v 
and a. It makes use of two independent four vectors, p and q, and the metric tensor. It contains 
fourteen Lorentz structures with fourteen invariant functions i^. Although all the Lorentz structures 
are independent, the fourteen invariant functions can not be independent due to current conservation. 
Enforcing current conservation reduces the number of independent invariant functions. Usually we may 
take advantage of current conservation and write sum rules which are simpler. The best example is 
the two-point function of the vector mesons p and J/ip where, instead of studying separately the two 
possible independent structures g^ v and q^q^, we may use a combination of them which is manifestly 
conserved: 

IV(g) = U(q 2 )(- 9lMU q 2 + q^), (65) 

since q^H^u = 0. The sum rule is thus written for the scalar function II(g 2 ). Inspite of technical 
advantages of using (65), we may still choose to work with one single structure, either g^ v or q^q u . In 
the case of (64) current conservation yiels identities such as 

p"IW(p,p') = 0, (66) 

or 

p' a r vai >M) = o, (67) 

depending on the momentum of the J ftp or p meson. With the help of expressions (66) or (67) we 
might rewrite the sum rules in terms of combinations of different structures. However there is nothing 
wrong in working with individual Lorentz structures, as done in [20] and [23]. 

In Fig. 13 we show the typical difference between results obtained with different structures. In 
the Figure we see the form factors of the pDD vertex, where we have only two structures p^ and p' . 
The dashed line shows the case of an off-shell p meson. In this case the two structures give the same 
result. The solid and dot-dashed lines refer to the case of an off-shell D meson and they show the same 
quantity computed in the two different structures. We can observe a clear difference between the two 
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Fi gure 9: Pole (solid line) and continuum (dashed line) contribuition to QpD*D* (Q 2 = lGeV 2 ), 
function of the Borel mass M 2 . 

lines, especially in the low Borel mass region. Therefore the choice of structures deserves attention, 
since it may be an additional source of uncertainties. In Fig. 14 we show the form factors of the 
J/ijjD*D* vertex for fixed M 2 as a function of Q 2 . The dashed and dotted lines show results obtained 
with different structures. 

4 Form factors 

In the previous section we have presented the conditions which must be satisfied for a sum rule to 
be considered reliable. In this section we discuss a difficulty inherent to the calculation of coupling 
constants with QCDSR. The solution of (31) is numerical and restricted to a singularity-free region in 
the Q 2 axis, usually located in the space-like region. Therefore, in order to reach the pole position, 
Q 2 = —m 2 Ms , we must fit the solution, finding a function g^M2M 3 (Q 2 ) which is then extrapolated to 
the pole, yielding the coupling constant. In the following subsections we introduce an extrapolation 
procedure and discuss how to improve it. 

4.1 The extrapolation procedure 

In order to minimize the uncertainties associated with the extrapolation procedure, for each vertex we 
perform the calculation twice, putting first one meson and then another meson off-shell, obtaining two 
form factors g^M 2 M 3 

(Q 2 ) and 

9^MiM 2 Mz(Q 2 ) an d requiring that these two functions have the same value 
at the respective poles. The superscripts in parenthesis indicate which meson is off-shell. 

In order to yield reliable results the sum rule (31) must satisfy the quality criteria discussed in the 
previous subsections. In first place (31) is a function of two Borel masses M 2 and M' 2 . A good sum 
rule is independent of the choice of these masses (it is "Borel stable"), showing a plateau when plotted 
as a function of M 2 or M' 2 . Moreover, the OPE side is a series which must be convergent. We choose a 
value for M 2 and plot (31) as a function of Q 2 , as it is shown in Fig. 15. The squares and circles show 
the result of the numerical evaluation of the form factor g(Q 2 ) as a function of Q 2 for the D*Dir vertex 
with a pion (squares) and a D (circles) off-shell. As it can be seen, at a certain (low) value of Q 2 the 
calculation stops, because at this point the stability and convergence criteria are no longer satisfied. 
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Figure 10: Pole versus continuum contributions to g p £)*D* 

(Q 2 = 1 GeV 2 ) function of the Borel mass 

M 2 . 

From here on we have to extrapolate. In Fig. 15 we show fits of the QCDSR results represented by the 
lines. As it is illustrated in the case of an off-shell pion, the numerical points obtained with QCDSR 
can be fitted, with similar accuracy, by several forms which, when extrapolated to Q 2 = —m 2 , will lead 
to very different points! This can be clearly seen by comparing the dashed and dash-dotted lines in 
Fig. 15. In order to reduce the freedom in the extrapolation and constrain the form factor we calculate 
and fit simultaneously the values of g(Q 2 ) of the same D*Dn vertex with the D off-shell. The results 
are shown in Fig. 15 with circles fitted by the solid line. We perform the fits of the two sets of points 
(circles and squares) imposing the condition that the two resulting parametrizations, when extrapolated 
to Q 2 = —m^ and Q 2 = —m 2 D go to the same value of gD*Dn(Q 2 )- This procedure is enough to reduce 
the uncertainties and, imposing this requirement leads to gD*Dn = 14.0 ± 1.5, which is consistent with 
the experimental value go*D-K = 17.9 ± 0.3 ± 1.9. 

As another interesting example, we consider the vertex J/ipDD*. Here we try to improve the 
procedure described above, calculating three form factors (one for each off-shell particle). This new 
procedure imposes a more stringent condition on the parametrizations. Fixing M 2 and M' 2 to the 
values of the incoming and outgoing meson masses we show, in Fig. 16, the momentum dependence 
of the QCDSR results for the three form factors g[j/^ DD *, gf/^DD* anc ^ #j/?/>dd* through the circles, 
squares and triangles respectively. Since our approach cannot be used at Q 2 <C 0, in order to extract 
the gj/^DD* coupling from the form factors we extrapolate the curves to the mass of the off-shell meson, 
shown as open circles in Fig. 16. In order to do this extrapolation we fit the QCDSR results with an 
analytical expression. We tried to fit our results to a monopole form, since this is very often used for 
form factors, but the fit was only good for gj J /^ D *- For g j% DD * and g^f/J DD * we obtained good fits 
using a Gaussian form. These fits are also shown in Fig. 16 through the dotted, solid and dashed lines 
respectively. From Fig. 16 we see that all three form factors lead to compatible values for the coupling 
constant when extrapolated to the off-shell meson masses (open circles in Fig. 16). 

From the parametrizations we can also extract the cutoff parameter, A, associated with the form 
factors. The general expression for the Gaussian parametrization is: Aexp [— (Q 2 + B) 2 /A 4 ], which 
gives A ~ 4.5 GeV when the off-shell meson is D or D* . For the monopole parametrization the general 
expression is: g[(A 2 — m 2 )/(A 2 + Q 2 )}. Therefore, for an off-shell J/ip we get A ~ 7.5 GeV. It is very 
interesting to notice that the value of the cutoff is directly associated with the mass of the off-shell 
meson in the vertex. The form factor is harder (i.e., the curve in the Figure is flatter) if the off-shell 
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Figure 11: Continuum threshold dependence of the form factor g!^j D * (Q 2 )- The dashed solid and dotted 
lines give the parametrization of the QCDSR results for A s = 0.4 GeV, 0.5 GeV and 0.6 GeV respec- 
tively. The lower, intermediate and upper set of curves show the results for A u = 0.4 GeV, 0.5 GeV 
and 0.6 GeV respectively. 

meson is heavier. 

4.2 Hadronic loops 

Coming back to the D*Dir vertex, we observe that, while the obtained number, go*Dn = 14.0 ± 1.5, is 
not far from the experimental value, there is still a discrepancy. The procedure of fitting the QCDSR 
points in the deep euclidean region and extrapolating them to the time-like region contains systematic 
uncertainties associated to the analytical form chosen for the parametrizations, i.e., monopole, expo- 
nential or gaussian. We tried to reduce this systematic uncertainty performing a double (and also a 
triple as in the case of the J/ijjD*D vertex) fit. However it would be desirable to have a more physical 
way to reduce the systematic uncertainty. 

In [37] the authors revisited this problem, employing hadronic loops, calculated by means of effective 
field theories (EFT), in order to produce a better parametrization for D*D7r results calculated with 
QCDSR. Purely hadronic calculations are independent from QCDSR and involve the choice of an 
effective Lagrangian, including the possible requirements of chiral symmetry and/or SU(4). Beyond 
the tree level, one has to deal with the problems and uncertainties associated with renormalization. 
As it was discussed in [37] a suitable combination of EFT and QCDSR results allows the reduction of 
undesired indeterminacies of both approaches, improving their predictive power. 

The full D*D7r vertex function in a hadronic approach involves the computation of several diagrams. 
Leading contributions to this vertex come from both the tree interaction and from the diagrams depicted 
in Fig. 17. Meson loops are a necessary consequence of quantum field theory and do contribute to several 
hadronic observables. In practice, due to problems associated with infinities, renormalization becomes 
unavoidable in the evaluation of loop corrections to observables. Here the basic idea is to isolate the 
unknown loop parameters into some basic constants, in such a way that they can be determined by 
matching the results of loop and QCDSR results. 

We note that some diagrams, such as, for instance, that in Fig. 17 a, contain internal vertices 
involving the D*Dn coupling. This suggests that the calculation is "cyclic", since one needs to use the 
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Figure 12: Dependence of the D*Dp form factors on the continuum thresholds. The steeper lines 
correspond to p off-shell. The others are for D off-shell. 

D*D7r form factor in order to calculate the D*D7r form factor. Actually, there are differences between 
the internal particles and the external ones. The former are always virtual, whereas the latter may 
be either real or put on mass shell in the extraction of the coupling constant. In the framework of 
perturbation theory, at leading order, internal particles are treated as elementary, without structure. 
They are assumed to be point-like and the evaluation of leading terms does not require the use of 
internal form factors. Consistently, one must use bare coupling constants for these interactions. 

Since there are heavy mesons circulating in the loops shown in Fig. 17, one might argue that other 
states should also be included. We do have, for example, fermion-antifermion components such as NN 
or A C A C in the loops. An incoming positive pion can split into a p plus a n, and so on. However, in a 
different context [38], it has been shown that this kind of splitting is suppressed with respect to the pion 
— > meson-meson splitting, by one order of magnitude. The neglect of this kind of contribution seems 
therefore justified. The same holds for the possibility of strangeness circulating in the loop, associated 
with virtual states such as D s , D*, K and K*. Using only 7r's, p's, D's and .D*'s the low and high 
Q 2 regions of the form factor are covered. Thus it is enough to work with a simple effective theory. 
It is convenient to use the effective Lagrangian (40), which is constrained by SU(2) flavor and chiral 
symmetries, as well as gauge invariance. The coupling constant appearing in the Lagrangian is the bare 
one. The p couplings are assumed to be universal and are implemented by covariant derivatives of the 
form V th = — ig p T p^ , where g p is the universal coupling constant and T is the isospin matrix 
suited to the field upon which the derivative acts. With the proper Lagrangians it is possible to 
write and evaluate all the contributions to the total vertex function. 

The 7i(q) D{p) D* a {p') vertex function T p (p 2 ) for an off-shell D is written as: 

r> 2 ) = - q.g^Lip 2 ), (68) 

where g^Dnip 2 ) is the form factor, such that the physical coupling constant is gD*Dw = fi , D*D 7 r( m i))- 
In [37] two kinds of loop corrections to this vertex were considered, containing pion and p intermediate 
states, denoted respectively by F n (p 2 ) and F p (p 2 ). The perturbative evaluation of these functions gives 
rise to divergent integrals and gD^DwiP 2 ) can be determined only up to yet unknown renormalization 
constants. The use of standard loop integration techniques, such as dimensional regularization and MS 
subtraction of divergences, for all diagrams, allows one to write the form factor as: 

</Sk(p 2 ) = K + C n F^p 2 ) + C p F p (p 2 ) , (69) 
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Figure 13: M' 2 dependence of the DDp form factors at Q 2 = 1 GeV 2 for A s = A u = 0.5 GeV. The 
dashed line gives the QCDSR result for g^DpiQ 2 ) an d the dot-dashed and solid lines give the QCDSR 
results for g^DpiQ 2 ) m the p^ and p'^ structures respectively. 



where K , C n and C p are constants. These constants incorporate the bare couplings, the usual parameters 
associated with renormalization and here they are determined by comparing g^D^ip 2 ) with the results 
from QCD sum rules. Keeping only the terms which depend on p 2 the explicit evaluation of the 
diagrams can be performed and the form factor as a function of p 2 , the D four-momentum squared, 
can be obtained. At this stage, it still contains three unknown parameters, which are determined 
by adjusting the function g^Diriv 2 ) to the QCD sum rule points taken from [19]. Those results are 
displayed in Fig. 18, where P 2 = —p 2 , together with the best fit (x 2 ~ 1CT 3 ) represented by the solid 
line. Computing the value of g^} D7T (p 2 ) at p 2 = m 2 D , one arrives at the following value for the coupling 
constant: 

g D * m = 17.5 ±1.5, (70) 

in very good agreement with experiment. The errors quoted come from the QCDSR points, which 
contain a typical error of ~ 10 %. In the same Figure we can also see the results of the fits of the 
QCDSR points with two mixed monople-dipole structures with three free parameters, namely: 



g\p 2 ) 



g n (p 2 



c 



Ci 



A 2 — m 




A 2 



m 



D 



A 2 



+ C 2 



'A 2 - m 2 D 



A 2 



P 



(71) 



(72) 



which yield xj ~ 10~ 3 (dashed line) and xji ~ 10~ 2 (dash-dotted line), respectively. Looking at Fig. 
18 we learn that these alternative structures, reasonable as they are, diverge significantly from the loop 
calculation in the region where the D is not too off-shell, stressing the importance of a proper hadronic 
treatment of the form factor in that region. 

As far as practical applications are concerned, our numerical results for the form factor gD*DTr{p 2 ), i n 



the whole range 



-m 



D 



< P 2 < 5 GeV 2 , are very well described by the mixed monopole-dipole structure 



given by Eq. (71) with the parameters C = 8.7, A x = 5.1 GeV and A 2 = 2.9 GeV. These results suggest 
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Figure 14: gj/^ D * D * (circles and squares) and gj/^* D * (triangles) form factors as a function of Q 2 
from the QCDRS calculation of this work. The solid and dashed lines correspond to the monopole 
parametrization of the QCDSR data, using two different structures, and the dotted line corresponds to 
the exponential parametrization. 

that the use of meson loops can reduce the uncertainty in the extrapolation of form factors, computed 
in the space-like region by means of QCDSR, to the time-like region, with the corresponding increase 
in the reliability of predictions for coupling constants. Apart from the approximations described above, 
the procedure has no new source of errors. 

To conclude this section, we have discussed a new method of improving QCDSR calculations of 
hadronic form factors, which consists in matching QCDSR results, valid mainly in the deep euclidean 
region, to meson loop calculations, valid when the D is not too off-shell. This matching is well justified 
from the physical point of view, since in the intermediate and large Q 2 regions the relevant degrees of 
freedom are the quarks and gluons, with non-perturbative corrections taken into account through the 
QCD condensates. The opposite happens for low values of Q 2 , where sum rules calculations become 
non- reliable due to the lack of a large mass scale. At this point, the meson exchange dynamics becomes 
a reliable tool, but it depends on unknown constants associated with the renormalization of the mesonic 
vertices. Although the exact frontier between meson dynamics and QCDSR cannot be precisely known, 
the success of the method in the example considered here supports the view that the matching may 
become useful in increasing the predictive power of both procedures. 



5 Results 

5.1 Form factors and couplings 

We have applied the procedure described in the previous sections to the following vertices D*D7r, 
D*D*n, DDp, D*Dp, D*D*p, J/ipDD, J/i)D*D and J/i)D*D*. As mentioned above, for each vertex 
we obtain two sets of points, which have been parametrized by the following forms: 

{I) (73) 
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Figure 15: D*Dir form factor as a function of Q 2 for an off-shell pion (dashed and dot-dashed lines) and 
an off-shell D (solid line). Circles and squares are the results of numerical calculations and the curves 
are their fits. 

(U) 9 { mXm 3 = Aex V [-(Q 2 /B)] (74) 
(HI) 9 ( m;Lm 3 = AeM-{Q 2 + Cf/B}. (75) 

In Table 3 each line refers to the vertex indicated in the first column. In the second, third and fourth 
columns we present the values of the parameters A and B for the case where a heavier meson in the 
vertex (Mi) is off-shell, indicating also which parametrization ((I), (II) or (III)) was employed. In the 
fifth, sixth, seventh and eighth columns we show the parameters and type of parametrization used in 
the case where a lighter meson (M 2 ) is off-shell. The numbers presented in Table 4 summarize our 



Mi M 2 M 3 






Mi off 






M 2 off 






Form 


A 


B 


Form 


A 


B 


C 


DnD* 




126 


11.9 


(II) 


15.5 


1.48 




J/ipDD* 




200 


57 


(HI) 


13 


450 


26 


J/ipDD 




306 


63 


(HI) 


15 


250 


20 


DpD 




37.5 


12.1 




2.5 


0.98 




J/ipD*D* 




400 


78.5 




1.96 


3.5 




D*pD* 


(II) 


4.9 


13.3 




5.2 


2.7 




D*nD* 


(II) 


4.8 


6.8 




8.5 


3.4 




DpD* 


(I) 


234 


44 




5.1 


4.3 





Table 3: Parameters used in (73), (74) and (75). B and C are in GeV 2 and A is either in GeV or GeV 2 , 
depending on the vertex. All the isovector mesons it and p appearing in the Table are charged. 

results. They contain uncertainties coming from different sources, which will be discussed in the next 
section. The vertex coupling constants can be obtained from the form factors and they are presented 
in Table 4. A comparison with other estimates is going to be done in section 6. 
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Figure 16: Momentum dependence of the J/ipDD* form factors. The dotted, dashed and solid lines 
give the parametrization of the QCDSR results (triangles, squares and circles). 



M 1 M 2 M 3 


9 Mi M 2 M 3 


DttD* 


9.9 ± 1.0 


J/ipDD* 


4.0 ±0.6 GeV- 1 


J/ipDD 


5.8 ±0.9 


DpD 


3.0 ±0.2 


J/ipD*D* 


6.2 ±0.9 


D*pD* 


4.7 ±0.2 


D*nD* 


6.1 ±0.7 GeV" 1 


DpD* 


4.3 ± 0.9 GeV" 1 



Table 4: Coupling constants for neutral isovector mesons (n° and p°). 



5.2 Uncertainties 

We consider now, as an example, the vertex D*Dp and discuss, one by one, all the sources of uncertainties 
in the calculation of the form factors. The sum rules for D, p and D* off-shell are given by: 

(D) 2 2 

n 9d*d p -^f 1 f , f , (d) ( + n —fr 
Cttt^ ^-e M' 2 e a/2 — / rf s / du p y >(u,s,t)e e m> 2 

(Q 2 + m 2 D ) Air 2 J J 



and 



-(qq)e-</ M \ (76) 

Q rr ^b m x)* 1 r r u 

C (Q2*+ P m 2 j e ~ 1777 e ~ = y rfs / du p {p) {u,s,t)e~^e~^, (77) 

( D *) m 2 m 2 J 

C (g'±mi) e "^ e "^ = -^ 2 \ ds \ du piD * ){u > s ' t} e "^ e "^ 

±(gg)e-^/ M ' 2 , (78) 
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Figure 17: Meson loop contributions to the D*Dir form factor: "triangle" diagrams. In the internal 
triangles the solid, dashed, wavy and spring lines represent a virtual D, n, p and D* respectively. 

where t = —Q 2 and the functions p^ D \ p^ and p( D *^ are the double discontinuities associated to the 
perturbative diagram of the sum rules with an off-shell D, p and D* off-shell respectively. They are 
given by: 



P {D \u,s,t) 



3m c 

71 



u(2m c — s — t + u) 



A 



(79) 



with A = (u + s — t) 2 — Aus. The integration limits in (76) are: 

m 2 (s + t — ml) — st 



< u < 



mi 



and 



m c < s < So, 

The perturbative contribution for the double discontinuity for an off-shell p meson is given by: 

3m r t 



A 3 / 2 



u + s — t — 2mi 



and the corresponding integration limits in (77) are: 



m 2 (s — t — ml) 
s — ml 



< U < Uq, 



and 



m c < s < s , 

The perturbative contribution for the double discontinuity for an off-shell D* meson is given by: 

3m r 



P 



(£>*). 



u, s, t) 



A 3 / 2 



s{2m 2 + s — t — u) 
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Figure 18: The D*Dn form factor. Dots: QCDSR from [19]; solid, dash and dash-dotted lines are fits 
obtained with eq. (69), (71) and (72), respectively. 

The integration limits in (78) are: 

m 2 {t — s — m 2 ) 



t < u < 



t — ml 
and 

< s < s , 

In the above expressions C is a constant defined as: 

n m r>fD , , 
C = m p f p m D *f D *, 

where fp, f P and fo* are the decay constants of the mesons D, p and D* respectively. The numerical 
evaluation of the above expressions is done with the numerical inputs shown in Table 2 and the resulting 
points are shown in Fig. 19. The triangles, squares and circles are the results for the g^DpiQ 2 )^ 

go*D P (Q 2 ) an d 9d*d p (Q 2 ) form factors respectively. As indicated in Table 3, in the case of an off-shell 
D meson, our numerical results can be fitted by the following monopolar parametrization (shown by 
the dashed line in Fig. 19): 

9M) = q^JJ. (82) 

where the function g^DpiQ 2 ) nas the units of GeV~\ Following the procedure discussed above, we de- 
fine the coupling constant as the value of the form factor at Q 2 = —m 2 M , where tum is the mass of the me- 
son M. Therefore, using Q 2 = —m 2 D in Eq (82), the resulting coupling constant is g^Dp = 5-76 GeV -1 . 
For an off-shell p meson our sum rule results can be fitted by an exponential parametrization, which is 
represented by the solid line in Fig. 19: 

^D P (Q 2 ) = 5.1e-« 2 /4.3. (g3) 

Using Q 2 = —m 2 in Eq (83) we get g^Dp = 5-89 GeV^ 1 . In the case of an off-shell D* meson, our 
numerical results can be fitted by the following monopolar parametrization (shown by the dotted line 
in Fig. 19): 

- TO' (84) 
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Figure 19: The D*Dp form factor. g D * Dp (squares), g£*Dp (triangles) and g D *^ p (circles) form factors 
as a function of Q 2 . The dotted, solid an dashed lines correspond to the parametrizations discussed in 
the text. The vertical bars show the theoretical errors in the coupling constants once all variations in 
the parameters are taken into account, as explained in the text. 



Evaluating this form factor at Q 2 = —m 2 D , we find the coupling g D *^ p = 6.65 GeV 1 . 

Looking at Fig. 19 we can observe that the D off-shell form fator is much harder than the p off-shell 
one. This agrees with the behavior observed in Figs. 14, 15 and 16: the heavier is the off-shell meson, 
the harder is its form factor. Following this same trend, we would expect the D* off-shell form factor to 
be even harder than the D off-shell one. However, comparing the dashed and dotted lines in Fig. 19, 
this seems not to be the case: the slope of the D* curve is slightly bigger than the one of the D curve. 
Since their mass difference is relatively small (~ 150 MeV) the two curves should have almost the same 
slope. The observed difference is an indication of the limited precision of our method. 

The form factors (82), (83), (84) and their extrapolations to the on-shell points leading to the 
coupling constants do not contain error bars. In fact, a careful and systematic study of errors in 
QCDSR calculations is hard to find in the literature. We took Refs. [39] as a guide. In Fig. 19 we can 
see the theoretical error bars at the endpoints of the three curves. In what follows we describe how we 
obtain them. First, we compute the sum rules (76), (77) and (78) extensively, taking into account the 
errors in the masses, decay constants, condensates, choice of the Borel mass and continuum threshold 
parameters. In each computation all the parameters are kept fixed, except one, which is changed 
according to its intrinsic error. The errors in the quark and gluon condensates, in the masses and decay 
constants are listed in Table 2. The three Borel masses were chosen in the interval 2.7 < M 2 < 3.3 
GeV 2 for an off-shell p and an off-shell D* and in the interval 27 < M 2 < 33 GeV 2 for an off-shell D. 
After each round of calculation of the three sum rules, we obtain three sets of points which are then 
fitted and extrapolated to the respective on-shell points. The sets of points are all fitted with the forms 
(82), (83) and (84), but for each set the numerical constants appearing in these forms are different. 

Every extrapolation introduces some ambiguity in the final results, since we have the freedom to fit 
a set of points with different parametrizations. In our case this freedom is strongly reduced because we 
require that all the three parametrizations lead to approximately the same coupling constant. In Fig. 
19 this requirement forces the three endpoints of (82), (83) and (84), which are taken at the squared 
masses of the corresponding particles, to coincide, i.e., to have approximately the same height in the 
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Quantity 




a 


a(%) 


A 


5.86 


0.08 


1.4 


fp 


5.89 


0.01 


0.1 


fo 


5.41 


0.65 


12.0 


Id* 


5.90 


0.40 


6.8 


M 2 


5.90 


0.10 


1.7 


m c 


5.97 


0.40 


7.4 



Table 5: Changes in g^l Dp induced by changes in different quantitities. 

figure. Of course, due to the approximations used, we can not expect this matching to be perfect. Once 
this procedure is completed and we determine the three coupling constants with an error corresponding 
solely to the variation of one parameter, we move to the next parameter to be varied, keeping all others 
fixed and repeat the procedure. In each step we can have an idea of how sensitive is each coupling 
constant to the parameter under consideration. In the end, for each coupling constant we take the 
average of all encountered values and calculate also the global error, which is shown in Fig. 19 as 
an error bar at the on-shell point. The final number is then obtained taking the average of the three 
couplings found and the final error is also obtained from the errors of each coupling. 

Among the sources of errors, one deserves a special discussion. Very often in QCDSR calculations, 
appreciable uncertainties in the results come from the lack of knowledge on the continuum threshold 
parameters. In order to study the dependence of our results on these parameters, we vary A SjU between 
0.4 GeV < A SjM < 0.6 GeV in the sum rule (77), between 0.4 GeV < A, < 0.6 GeV and 0.65 GeV < A u < 
0.75 GeV in the sum rule (76) and between 0.65 GeV < A s < 0.75 GeV and 0.50 GeV < A u < 0.70 GeV 
in the sum rule (78). This variation produces new sets of curves which are shown in Fig. 12 and give us 
an uncertainty range in the resulting coupling constants g^*Dp an d 9d*Do- F° r the sa ^e of clarity we did 
not include the lines corresponding to the coupling g^*D P - Surprisingly, in the case of the form factor 
9d*Doi we observe a convergence of the extrapolation lines, which reduces the final error. Due to this 
accident, the continuum threshold parameters are not, in the gD*Dp case, the ultimate source of error. 
Their contribution (denoted by A in the Tables) is still significant, as it can be seen in Tables 5, 6 and 
7, but now they have less impact on the final error than the uncertainties in the decay constants fo 
and fo* and in the charm quark mass. In these Tables we show in the first column the quantity which 
was varied, in the second the average coupling constant resulting from that variation, in the third the 
standard deviation and in the fourth the percentual significance of a. 

After scanning the space of reasonable values of all the parameters, we conclude that, in spite of 
the inherent uncertainties, the sum rules really point to a value of the coupling constant! Of course, 
as in most of QCDSR calculations, the lack of precision is due to the "usual suspects", i.e., continuum 
threshold parameters, decay constants, heavy quark masses and condensates. A comparison of the 
Tables shows an intriguing aspect, namely that some of the input quantities affect each of the three 
sum rules in a quite different way. This may be a signal that some of the sum rules are less reliable than 
others. A deeper investigation of this question would involve several refinements, such as the calculation 
of a s corrections and higher order terms in the OPE. 

Considering the results presented in the Tables, the couplings are: 

g W Dp = 5.71 ±0.62 GeV" 1 , 
g $l Dp = 5.87 ±0.53 GeV" 1 , 

and 

ygf) = 6.63 ± 0.73 GeV" 1 , 
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Quantity (g [ pJ Dp ) a a (%) 



A 


5.95 


0.87 


14.7 


fn 
J p 


5.76 


0.01 


0.1 


fn 
J u 


5.30 


0.64 


12.0 


fo* 


5.80 


0.40 


6.8 


M 2 


5.76 


0.05 


0.8 


m c 


5.70 


0.30 


5.6 


(qq) 


5.77 


0.04 


0.8 



Table 6: Changes in g D J Dp induced by changes in different quantitities. 



Quantity 


a (D ) ) 

\yD*Dp/ 


a 


a(%) 


A 


7.00 


1.00 


14.3 


n 


6.61 


0.07 


1.1 


fn 


6.11 


0.74 


12.0 


I'd* 


6.69 


0.46 


6.8 


M 2 


6.65 


0.19 


2.8 


m c 


6.61 


0.06 


0.8 


(qq) 


6.66 


0.08 


1.3 



Table 7: Changes in g D ^ p induced by changes in different quantitities. 



We can see that the three cases considered here, off-shell D, p and D*, give compatible results for the 
coupling constant. Considering all the uncertainties and taking the average between the obtained values 
we have: 

9d*d p = (6.1 ±1.3) GeV-\ (85) 

Our results were obtained for certain concrete choices of currents, which represent charged states. 
Consequently, the obtained couplings are for charged states. As it will be discussed in the next section, 
from the coupling of charged states we can get the coupling of neutral p states, which will be identified 
with the "generic coupling", through the relation: 

9p+D°D*+ 9p-D»D*+ / oc x 
9D*Dp - 9D*Dp° - ^= - ^= • (°t)J 

Therefore the final value of the coupling constant listed in Table 4 is: 

g D » Dp = (6.1±1.3)/v^= (4.3 ±0.9) GeV -1 , 

To close this section we emphasize that, from the analysis of the Tables we conclude that, in the present 
context, the average error of our calculations is in the range from 10 to 15 %. 



6 Discussion 

In this section we compare our results with other QCDSR calculations, in particular with light-cone 
QCD sum rules (LCSR) results, and also with results obtained with other techniques. Besides QCDSR, 
coupling constants can be estimated with the vector meson dominance (VMD) model, with heavy quark 
effective theory (HQET), with SU(4) symmetry relations, with effective models such as the constituent 
quark- meson model, with chiral models and with lattice QCD calculations (LQCD). 
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6.1 SU(4) and HQET 



If SU(4) would be exact, several relations between the coupling constants should hold. Using the QCDSR 
results reported in Table 4 we can check to what extent these relations are satisfied. The relevant 
relations and their deviation from our results are shown in Table 8. They are ordered by increasing 
degree of violation. Remembering the typical uncertainties of 10 % to 15 % in our calculations, we can 
conclude that the first seven SU(4) relations are reasonably satisfied whereas the last six relations are 
badly violated. Although there is no rigorous systematics, we can clearly observe that violations occur 
mostly when there is a pion in the vertex. 



SU(4) Relation 


Violation 


g.j/i>DD 


— 9j/4>d*d* 


(7%) 


gpDD* = 


—gj/ipDD* 


(12%) 


9pDD = 


Ve>„ 

—g.j/i>DD 


(17%) 


9nD*D* - 


- —9j/ipdd* 


(20%) 


9D*D*p = 


— gj/i)D*D* 


(20%) 


9DDp = 


— g.J/i>D*D* 


(21%) 


9pD*D* - 


- —g.J/i>DD 


(25%) 


9nD*D 


' — g P DD* 


(29%) 


9pDD 


= gpD*D* 


(36%) 


9D*Dtv 


= gu*D*p 


(52%) 


9D*Dn = 


— g.j/i>D*D* 


(62%) 


9D*Dn = 


- —g.j/4>DD 


(64%) 


9d*Dtv — 9DDp 


(70%) 



Table 8: SU(4) relations between the coupling constants (on the left column) and their violation (in 
percentage on the right column) found in QCDSR. 



We can also check the heavy quark spin symmetry relations [40], which are presented below with 
their deviation from our results (calculated with the values presented in Table 4 in parenthesis: 

gpD*D = 9j ^- (45%) (87) 
m£,* 

gj/^D = (22%) (88) 

m D 

These relations come from HQET. Since charm is not heavy enough to ensure the validity of HQET, 
we would expect a significant violation of the above relations. In fact, considering the error in our 
calculations (88) is still satisfied whereas (87) is severely violated. 



6.2 Light cone sum rules 

In Table 9 we present a compilation of the estimates of the coupling constants go'Dn and gB*Bw from 
distinct calculations. From this Table we see that our result in Ref. [19] is in a fair agreement with 
the LCSR calculation in Refs. [42, 44], but is still smaller than the experimental value [49]: gwo-n = 
17.9 ±0.3 ± 1.9. However, using a better way to extrapolate the QCDSR results based on a meson loop 
calculation, the result obtained in Ref. [37] is in an excellent agreement with the experimental value. 
From this Table we also see that LQCD results are in a very good agreement with the experimental 
value. 
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Approach 


9D*Dit 


9B* Bit 


QCDSR [41] 


9±2 


20 ±4 


QCDSR [41] 


7±2 


15 ±4 


LCSR [42] 


11 ±2 


28 ±6 


QCDSR [43] 


6.3 ± 1.9 


14 ±4 


LCSR [44] 


10.5 ±3 


22 ±9 


QCDSR [19] 


14.0 ± 1.5 


42.5 ±2.6 


QCDSR plus meson loops [37] 


17.5 ± 1.5 


44.7 ± 1.0 


LQCD [45] 


20 ±2 




LQCD [46] 


18.811* 




dispersive quark model [47] 


18 ±3 


32 ±5 


Dyson-Schwinger equations [48] 


i5.8±?:S 


30.0±U 



Table 9: Summary of estimates for gD*Dw and gB*Bn- These couplings refer to charged mesons tt^. 



The basic difference between the QCDSR, described here, and the LCSR is the fact that, instead 
of considering the three-point function in Eq. (15), the central object in the LCSR is the correlation 
function of two meson currents between the vacuum and one on-shell meson state [50]: 

T( q ,p) = J d 4 x e* x (M 1 (p)|TO 3 (x)i 2 t (0)}|0). (89) 

The idea is to expand the product of the currents near the light-cone x 2 = 0. This expansion is different 
from the local OPE expansion used in the QCDSR because it incorporates a summation of an infinite 
series of local operators [50]. In the case that M\ is the pion, the vacuum-pion matrix elements are 
expressed via pion light-cone distribution amplitudes (DA). Since the pion DA's have a well defined 
twist, the LCSR can be written in terms of a twist expansion. On the other hand, if M\ is not light, 
its DA does not have a well defined twist. If M\ is a heavy meson, like a B meson, the correlation 
function can be systematically expanded in the limit of large m?, in heavy quark effective theory [51]. 
However, for hadronic vertices involving only charmed mesons there is no well defined expansion for the 
correlation function in the LCSR approach. An attempt to use LCSR to calculate the D*D*p coupling 
constant was performed in [52]. The obtained value is quoted in Table 10. It is closer to SU(4) and 
to VMD (see below) than the value found in [23]. However, a discussion of some potential sources of 
uncertainties, such as the use of a relatively large value of the Borel mass and the possible continuum 
dominance, is still missing in [52]. 



6.3 Vector meson dominance 

In the VMD model [2] a virtual photon with four-momentum q is emitted in the process eM — > eM, 
where M represents any meson. The photon can be decomposed into a sum of all neutral vector mesons 
including both isospin and isospin 1. Then the vector meson couples with the external meson M. At 
q 2 = one can write: 

£ 1X9VMM = ^ (go) 

v=p,u,,<M/W-. mv 

where •jy is the photon-vector-meson coupling that can be determined from the vector meson partial 
decay width to e + e~: 

r - ai2y ran 

^ " 3mf/ (91) 

Using the VMD model the couplings g P DD, 9pD*D*, Qj/ipDD and gj/^n*D* were estimated in Ref.[32]. 
The couplings g p o*D and gj/^D*D were estimated in Ref. [33] by applying the VMD model to the 
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Coupling 


QCDSR 


VMD 


Other models 


9pDD 


3.0 ±0.2 [17] 


2.52 [32] 




9pD*D (GeV- 1 ) 


4.3 ± 0.9 [24] 


2.82 [33] 


4.17 ± 1.04 [54] 


9pD*D* 


4.7 ±0.2 [23] 


2.52 [32] 


1.8 ±0.5 [52] 


QwDD 


-2.9 [53] 


-2.84 [32] 




9j/iPDD 


5.8 ± 0.9 [22] 


7.64 [32] 


8.0 ± 0.5 [40] 


9j/4,d*d (GeV- 1 ) 


4.0 ± 0.6 [22] 


8.0 ±0.6 [33] 


4.05 ± 0.25 [40] 


gj/i/jD*D* 


6.2 ±0.9 [20] 


7.64 [32] 


8.0 ± 0.5 [40] 



Table 10: Summary of estimates for gvD(*)D(*)- I* 1 the above couplings p stands for p . 

radiative decay D* — > Dj. The obtained values are presented in the Table 10 where we also present a 
summary of the predictions for the coupling constants in the Z)W vertex. 

In our approach the coupling constant is given by the value of the form factor at Q 2 = —m 2 M , where 
M is the off-shell vector meson. In the case of J/ip this corresponds to Q 2 = —9.6 GeV 2 whereas for 
Po the on-shell point is Q 2 = —0.6 GeV 2 . In the VMD model the vector meson has Q 2 = GeV 2 . 
We would then expect our results to present a reasonable agreement with the VMD estimates of the p 
couplings and a significant disagreement with the VMD estimates of the J/tp couplings. This pattern 
is not so evident in Table 10, which shows the approximate expected behavior in some cases but also 
shows discrepant behavior in other cases. All the form factors decrease as Q 2 goes from the time-like to 
the space-like region and hence all the values shown in the left column of Table 10 will become smaller 
when calculated at Q 2 = GeV 2 . This will improve the agreement with the VMD estimates in the 
case of the light vector mesons and will increase the disagreement with the heavy vector mesons. Since 
QCDSR are more reliable in the latter case, this comparison suggests that the use of VMD for heavy 
mesons is dangerous. 

7 J/ip absorption and production 

As an application of the form factors obtained above we address now the problem of J/ip absorption and 
production in hadronic matter. With the Lagrangians (40) - (50) we are able to compute the process 
DD — > J/ip + n, which involves the diagrams in Fig. 20, the process D*D — > J/ip + n, corresponding to 
the diagrams shown in Fig. 21 and also the process D*D* — > J/ip + 7r, corresponding to the diagrams 
in Fig. 22. 

As extensively discussed in previous works, although the above Lagrangians and amplitudes are 
quite satisfactory from the point of view of symmetry requirements, their straightforward application 
to the computation of cross sections leads to unacceptably large results. This comes from the fact that 
the exchanged particles may be far off-shell and therefore they enter (or leave) a vertex with a very 
different resolving power. In one extreme case, a virtual J/ip probing a D meson, may behave like a 
parton. Of course, when this happens, the compact J / ip almost misses the large D and as a consequence 
the cross section of the whole process drops significantly. This physics of spatial extension and resolving 
power is contained in the form factors. It has been realized by many authors that calculations with and 
without form factors lead to results differing by up to two orders of magnitude! Therefore we simply 
can not ignore the form factors. We must include them in order to obtain reliable results! 

Looking at the diagrams in Figs. 20, 21 and 22 we notice that we need the following form factors 
(and the corresponding coupling constants): gi D Dr \*{t), 9j}l DD (t), #j/Jdd*(*)> #j/2?*d*(*) 

and g nD J D *(t), where t is the usual momentum transfer squared and in the superscript in parenthesis 
we denote the off-shell particle. This is an important distinction, because the form factors in the same 
vertex, as we have seen, are very different when different particles are off-shell. 
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Figure 20: Diagrams which contribute to the process DD — > J/ip + ir. 



The cross sections for secondary J/ip production is related to the annihilation through detailed 
balance. In Fig. 23 we show the J / ip secondary production cross section as a function of a/s, without 
form factors. In all Figures, the channels DD — > J/ip + n, DD* — > J/ip + 7r and D*D* — >• J/ip + tt are 
represented by solid, dashed and dotted lines respectively. In Fig. 24 we show the corresponding inverse 
reactions. As it can be seen, the cross sections have the same order of magnitude in both directions. 
Figs. 25 and 26 are the analogues of 23 and 24 when we include the form factors in the calculations. Of 
course, only these last two Figures correspond to realistic numbers. The comparison of the two sets of 
Figures is interesting to estimate the effect of form factors. In previous studies doing the same kind of 
comparison, as for example in [33] , the introduction of form factors reduced the cross sections by factors 
ranging between 20 and 50 depending on the channel. In that work the form factor was the same for 
all vertices and the cut-off, not known, was estimated to be between 1 and 2 GeV. Our study is much 
more detailed and not only each vertex has its own form factor, but, depending on which particle is 
off-shell the form factor is different. The final effect of all these peculiarities is the reduction of the cross 
sections by a factor around 7. Although significant, this reduction is smaller than previously expected. 

Fig. 25 contains our main results. The plotted cross sections can be compared with the results of 
[55] and, more directly, with [56]. In Fig. 2 of [55], although the variables in the plot are different, we 
can observe the same trend and relative importance of the three channels. In that work, the results 
were obtained with the quark model of [57]. Our curves share some features with the results of [56], 
such as, for example, the dominance of the DD* channel and the falling trend of the DD* and D*D* 
channels. The behavior of the DD channel is quite different. In the energy range of y/s > 4.5 GeV 
our cross sections are smaller by a factor of 2 (DD*) or 5 (D*D* and DD). These discrepancies are 
large but they are expected since in [56] all channels include the final state J/ip + p, which we did not 
include. In the model used by the Giessen group [56] the cross sections for D + D — > J/ip + it and 
D + D — > J/ip + p are similar and the same conclusion holds for the other inital state open charm 
mesons. If this would remain true in the effective Lagrangian approach, then our results including both 
final states would come closer to those of [56], giving thus a more theoretical support to the model 
considered there. 

The exercise presented in this section was meant to illustrate the use of the charm form factors 
discussed in this review. The ultimate computation of charmonium interactions in a hot and dense 
medium should include other effects, not mentioned here. For recent works on the subject see, for 
example, [58, 59]. 
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D" D D" D D" D D" D 

(2a) (2b) (2c) (2d) 



Figure 21: Diagrams which contribute to the process D*D — > J/ip + ir. 




D" o D' D D' D D" o D" 

(3a) (36) (3<: ) (id) (3c) 



Figure 22: Diagrams which contribute to the process D*D* — > J/ip + ir. 



' 1 1 1 1 1 1 1 ' 1 1 1 1 1 


J~s{GeV) 






DD -> \|/+ti 




DD* -> V|/+7T 


i.i.i, 


D*D* ->u+~ 

i . i , i , i 



4~s(GeV) 

Figure 23: J/ip secondary production cross section without form factors. 
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100 




4~s(GeV) 



Figure 24: J/if> absorption cross section obtained through detailed balance without form factors. 



a(mb) 





I 1 I 1 


' 1 1 1 1 1 1 1 






DD -> \\i+n 






DD* -> v|/+ti 




I.I. 


D*D* -> V|/+7t 

i i i i 



4~s{GeV) 

Figure 25: J /if) secondary production cross section with form factors. 
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10 



■ (mb) 



■y+TC-> DD 
\\i+k-> DD* 
V(/+7t-> D*D* 



0.1 




4.4 4.6 4.8 

■Js(GeV) 



Figure 26: J /if) absorption cross section obtained through detailed balance with form factors. 

8 Summary 

We have studied the form factors of vertices with charm mesons. They are relevant to understand data 
from heavy ion collisions at RHIC and LHC, from B decays at BELLE and BABAR and, in the future, 
charm production at PANDA. We have described how to calculate them with QCDSR and presented 
the results of our calculations. The comparison of the obtained coupling constants with those obtained 
with other methods shows that all the numbers have the same order of magnitude and the discrepancies 
between them go from a few percent to a factor two. In every approach there are improvements to be 
made and more accurate results are expected in the future. There are still vertices which have not been 
studied, such as those with i] c . With the already available charm form factors it is possible to address a 
number of problems of phenomenological interest. One of them, namely, the production and absorption 
of J /if} in nuclear matter was discussed here and the huge effect of form factors was shown in detail. 
Moreover some of these form factors have been used in [13] to study resonances formed through meson 
exchange interactions. 

Most of the calculations follow standard procedures in QCDSR but the extrapolation techniques to 
compute the coupling constant were developed by our group. We calculated two form factors for two 
off-shell mesons in the considered vertex. The simultaneous extrapolation of these form factors allowed 
the determination of the coupling constant. While this requirement proved to be crucial for finding 
the coupling, the inclusion of a third form factor, obtained by putting the third meson off the mass 
shell, with the subsequent triple extrapolation to obtain the coupling, did not bring any significant 
improvement in the results. 

The extrapolation method used in our works has a systematic error which comes from the choice 
of the analytic form of the extrapolating functions. We considered only monopole, exponential and 
gaussian parametrizations. However there is no physical reason for choosing these forms. Using the 
D*Dir vertex for an exploratory study, we calculated the form factor computing the relevant hadronic 
loops. In this way there is no need to guess a particular form for the form factor. Since the hadronic loop 
calculation contains non-perturbative physics it should be reliable in low Q 2 domain. The introduction 
of hadronic loops improved a lot our quantitative analysis of the D*D7r vertex. However we think that 
it would be premature to include it as an obligatory complement to the QCDSR formalism. Some loop 
diagrams were neglected because they were assumed to be less important but this must still be proven. 
Moreover the effective Lagrangian theories for other mesons and for baryons are less well known. 

The disagreement between some of our results and the SU(4) predictions has still to be clarified. 



40 



The same is true for the estimates made with the help of light cone sum rules, heavy quark symmetry 
and vector meson dominance. From the perspective of QCDSR there is still room for improvements, 
as, for example, the inclusion of a s corrections and higher order condensates. An important extension 
of our program will be to systematically study vertices with charm and strange mesons. Some of them 
were already considered in Refs. [60, 61, 62, 63, 64]. 

Beyond individual technical improvements, we believe that the QCDSR community should make a 
joint effort dedicated to systematically compare results and methods and arrive at a consensus on the 
present status of the calculation of coupling constants. With this review we wish to take a step in this 
direction. 
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